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(54) Precede de detection conjointe 

(57) La presente invention concerne un procede de 
detection conjointe destine a etre mis en oeuvre dans 
un recepteur d'un systeme de transmission, tel qu'un 
systeme de transmission a acces multiple a repartition 
decodes. Ladite detection consiste en une etape decal- 
cul des signaux qui ont ete emis par lesdits emetteurs 
dudit systeme mis sous la iorme d'un vecteur x et ce ; k 
partir de I'ensembie des signaux recus des emetteurs 
dudit systeme mis sous la forme d'un vecteur y, de I'en- 
sembie des caracteristiques dudit systeme mis sous for- 
me d'une matrice A et eventueilement de ('ensemble 
des caracteristiques de bruit mis sous forme d'une ma- 
trice n, ladite etape de calcul consistent en une etape 
diversion d'une matrice Mfonction de la matrice sys- 
teme A une etape d'appiication de I'ensembie des si- 
gnaux recus y sur une matrice resultant du produit de 
la matrice inverse Af 1 par le transpose conjugue de la 
matrice A ladite etape d'inversion de la matrice M con- 
sistant en sa decomposition de Cholesky en une matrice 
triangulaire U ou L 

II est caracterise en ce qu'il consiste a determiner 
ladite matrice triangulaire de maniere iterative, chaque 
iteration consistant d'abord au calcul d'un ensemble de 
blocs appartenant a une ligne de blocs de ladite matrice 
puis, en tonction de la vaieur prise par un critere prede- 
termine caracterisant revolution du processus de calcul 
de ladite matrice decomposed comparee k une vaieur 
seuil, soit au passage a Iteration suivante soit a la du- 
plication de la ligne courante sur toutes les autres lignes 
restantes avec decalage a droite d'un bloc et, Eventuel- 
lement, tronquage des blocs en excedant. 
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Description 

[0001] La presente invention concerne un procede de detection conjointe destine a etre mis en oeuvre dans un 
recepteur d'un systeme de transmission, tel qu'un systeme de transmission du type a repartition de codes CDMA (Code 
Division Multiple Access : acces multiple & repartition de codes). 

[0002] De nombreux systemes de telecommunication en cours de developpement mettent en oeuvre les techniques 
CDMA. Ces techniques utilisent le partage en code des ressources frequencies, c'est-a-dire que plusieurs utilisateurs 
utilisent la meme bande de frequence en meme temps et sont differences par les codes d'etalement differents qu'iis 
utilisent. Un utillsateur distingue en general le signal qui lui est adress<§ dans le signal total regu en comparant celui ci 
avec son code d'etalement. 

[0003] Cependant la qualite du signal regu est diminuee par ('interference residuelle venant des autres utilisateurs 
appeiee MAI (Multiple Access Interference - Interference d'acces multiple). Divers techniques dont la detection con- 
jointe. plus particulferement utiiisee dans le cas de codes d'etalement synchrones : visent k utiliser la connaissance 
des codes d'etalement des autres utilisateurs afin de reduire cette interference residuelle des autres utilisateurs. 
[0004] On a represents a la Fig. 1 un modele de la liaison montante d'un systeme de transmission du type CDMA 
qui utilise un procede de detection conjointe. Ce systeme est prevu pour recevoir K utilisateurs pourvus chacun d'un 
emetteur 10 k (k ~ I a K). Chaque emetteur 1 0^ est prevu pour transmettre les donnees d* d'un utilisateur k par salves 
de N symboles de donnees, ces donnees etant alors notees d*^ \x\>x%) ou x£ represente aiors le n'^symbole emis 
par J'utilisateur k, Dans chaque emetteur 10*, les donnees cf< d'un utilisateur k sont soumises a un code d'etalement 
c* de longueur O bits alloue a i'utilisateur /cconsidere avant d'etre emises sur un canal 20^. Chaque canal 20* est 
caracterise par une reponse impulsionnelle variant dans le temps note© Met une largeurde bande frequentielie notee 
B c^ 1 ST* Les donnees sont alors transmises sur un canal 20 k . Pour ^estimation de la reponse impulsionnelle h* de 
chaque canal 20*. sont prts en compte L echantillons espaces les uns des autres d'un temps qui est egal a T c . 
[0005] Klein and all ont pu montrer que le signal y regu par le recepteur peut etre represente par un vecteur y de 
dimension (N x Q + L-1) dont les composantes sont representatives des signaux emis par chacun des K utilisateurs, 
eventueilement additionne a un bruit gaussien represente par un vecteur n egalement de dimension (NxQ+L-1). 
On peut done echre : 

Ax + n 

ou, 

est le signal emis avant etalement, 

est le signal regu par le recepteur 

est le bruit blanc gaussien additif, et 

A est une matrice de dimensions (N.Q + L-1) x (K.N) 
[0006] Chaque element de la matrice A, note A f] , est determine en fonction de ('estimation de la reponse impulsion- 
nelle de chacun des canaux d'ordre k. 

A t =[A l ' ( / J *,. ) /i*javec* = l,...,/: 
et en fonction de chacun des codes d'etalement d'ordre k : 
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c*=t* f c*,... f c*javec* = l,... > A' 



s [0007] Si I'on pose 

/- Q.(n-1) \~ w 

10 et 

j = k + K(n-1) 

is [0008] Un systeme de transmission du type CDMA peut etre modelise d'une maniere ou chaque element A^peut 
s'ecrire sous la forme : 

A, j * pour \..K t n~ 1...Net *v= 1 ..CM.-1 

A/j-s 0 sinon 

avec 

b" =c J •** =[6 l i ,6 2 %... J ^ fX . I ]pour* = l..i: 



[0009] Afin d'estimer les donnees transmises par chaque utilisateur k, on peut proceder, dans un recepteur 30, a 
une detection conjointe telle que celle qui est decrite par AXIein, G.K.Kaleh et P.W.Baier dans un article intitule "Equa- 
lisers for multi-user detection in code division multiple access mobile radio systems" paru dans IEEE de 1994. 
[001 0] Ces procedes de detection conjointe connus consistent essentieliement a resoudre i'equation matricielle don- 
nee ci-dessus en determinant le vecteur des symboles x en connaissant le vecteur des signaux detectes y la matrice 
30 systeme A et la matrice de bruit n. 

[001 1 ] On notera que la connaissance des matrices y A et n est faite par ailleurs par des methodes d'estimation qui 
sont connues de I'homme du metier et qui ne sont par consequent pas decrites plus en detail dans la presente des- 
cription. 

[0012] On notera encore que, notamment du fait que la reponse impulsionnelie de chacun des canaux 20^ varie 
35 avec le temps, c'est a chaque instant que les elements de la matrice A change de valeurs. 

[0013] La determination du vecteur des symboles x, par la methode exposee ci-dessus, utilise generalement un des 
deux criteres suivants qui sont, d'une part, le forcage a zero (zero forcing ou ZF) et, d'autre part, Perreur quadratique 
moyenne minimum (MMSE = Minimum Mean Square Error). Ces deux criteres menent respective ment aux solutions 
suivantes : 

40 x * (A+A)' 1 A* y pour le critere de forcage a zero 

ei 

x - (A+A+(a 2 f)riA+y pour le critere d'erreur quadratique moyenne minimum MMSE (I est la matrice unite, o 2 etant 
la variance des elements de la matrice de bruit n). 

[0014] Dans la suite de la description, I'exposant + applique a une matrice signifie qu'il s'agit de son transpose 
^5 conjugue. 

[0015] Ces equations peuvent egalement s'ecrire : 

x = M*Vy 

50 

[0016] Comme on peut le constater, la resolution de cette equation necessite done I'inversron d'une matrice, notee 
M, qui est soit egale au produit A+.A de la matrice systeme transpose conjugue A+ par la matrice systeme A, soit egale 
au produit A+.A+<5 2 I de la matrice systeme transpose conjugue A+ par la matrice systeme A auquel est ajoutee la 
matrice unite I multiple par la variance o 2 du bruit. 
55 [0017] Lorsqu'eile est hermitienne, definie et positive, (ce qui est toujours le cas avec la moderation d'un systeme 
COMA donnee ci-dessus) une matrice Mpeut §tre decompose, au moyen d'un algorithme, dit algorithme de Cholesky, 
en une matrice triangulaire superieure U(ou inferieure L selon la notation choisie) telle que la matrice /West egale au 
produit de la matrice triangulaire superieure transposee conjuguee U*- par la matrice U : 



3 



EP 1 098 464 A1 

M= U*. U 

ou selon le cas, au produit de la matrice triangulaire inferieure L par fa matrice transposee conjuguee L + : 

5 

[0018] Or, I'inversion de la matrice M est grandement simplifiee puisque la matrice inverse At ? est soil egale au 
10 produit de la matrice inverse U' r par le transpose conjugue de cette matrice (If 1 p : 

M' 1 =U 1 (U 1 ) 4 

15 soit au produit de la matrice inverse conjuguee (L 1 )+ par la matrice L~ 1 : 

M' 1 ~(L* 1 f L' 1 

2Q [0019] Du fait du peu d'operations necessaires k inversion de U, Inversion de M en passant par la decomposition 
de Cholesky necessite beaucoup moins d'operations que I'inversion de M par un calcul direct 
[0020] Cette methode par decomposition de Cholesky autorise une reduction de la complexity de ce calcul par un 
f acteur 2 environ (complexity en T^/S, T etant la dimension de la matrice) par rapport aux techniques ^inversion directe 
classiques. Cependant la quantite de calcul reste encore tres importante. 
25 [0021] Le but de ia presente invention est de proposer un precede de detection conjointe qui utilise la methode par 
decomposition de Cholesky mais dont la quantite de calcul est dirninuee par rapport aux precedes de Tart anterieur 
[0022] A ce! effet, un procede de detection conjointe selon ia presente invention est du type qui consiste en une 
etape de calcul des signaux qui ont ete emis par lesdits emetteurs dudit systeme mis sous la forme d'un vecteur xet 
ce, a partir de Pensemble des signaux regus des emetteurs dudit systeme mis sous la forme d'un vecteur y, de Pen- 
so semble des caracteristiques dudit systeme mis sous forme d'une matrice A et eventueflement de Pensemble des ca- 
racteristiques de bruit mis sous forme d'une matrice n t ladite etape de calcul consistant en une etape ^inversion d'une 
matrice Mfonction de la matrice systeme A, une etape d'application de Pensemble des signaux regus ysur une matrice 
resultant du produit de la matrice inverse M' 1 par le transpose conjugue de la matrice A, ladite etape ^inversion de la 
matrice M consistant en sa decomposition de Cholesky en une matrice triangulaire Uou L 
35 [0023] Ce procede est caracterise en ce qu'i! consiste a determiner ladite matrice triangulaire de maniere iterative, 
chaque iteration consistant d'abord au calcul d'un ensemble de blocs appartenant a une ligne de blocs de ladite matrice 
puis, en fonction de la valeur prise par un critere predetermine caracterisant devolution du processus de calcul de 
ladite matrice decomposee comparee a une valeur seuil, soit au passage a Iteration suivante soit a la duplication de 
la ligne courante sur toutes les autres lignes restantes avec decalage h droite d'un bloc et ? eventuellement, tronquage 
40 des blocs en excedanl. 

[0024] Selon une autre caracteristique de Pinvention, ladite valeur seuil est determinee dynamiquemenl en fonction 
des ressources disponibles au moment du calcul. 

[0025] Selon une autre caracteristique de ['invention, ledit critere est tel que sa valeur depend des valeurs respec- 
tivement prises par au moins un element de Pensemble de blocs de la ligne courante et le ou les elements correspon- 
ds dants de Pensemble de blocs de fa ligne calculee & Piteration precedente. 

[0026] Selon une autre caracteristique de Pinvention, ledit critere est tel que sa valeur est egale a la valeur absolue 
de la difference des valeurs prises respectivement par I'element de Pensemble de blocs de la ligne courante dont la 
valeur absolue est la plus grande et par ('element correspondant de Pensemble de btocs de la ligne precedente. 
[0027] Selon une autre caracteristique de I'invention, ledit critere est tel que sa valeur est egale a la valeur maximaie 
50 des valeurs absolues des differences entre les valeurs prises respectivement partous les elements de I'ensembie de 
blocs de la ligne courante et par les elements correspondants de I'ensembie de blocs de la ligne precedente. 
[0028] Selon une autre caracteristique de i'invention, ledit critere a pour valeur le nombre d'iterations, ia valeur seuil 
correspondant aiors a un nombre maximal d'iterations. 

[0029] La presente invention concerne egalement un recepteur qui est prevu pour mettre en oeuvre un procede de 
55 detection tel qu'i! vient d'etre decrit. 

[0030] Les caracteristiques de ('invention mentionnees ci-dessus, ainsi que d'autres, apparaitront plus clairement & 
la lecture de la description suivante d'exemples de realisation, ladite description etant faite en relation avec les dessins 
joints, parmi lesquels : 
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La Fig. 1 est un schema synoptique d'un systfcme CDMA qui est prevu pour mettre en oeuvre un procede de 
detection conjointe selon la presente invention, 

La Fig. 2 est un schema d'une matrice resultant d'une decomposition de Cholesky mettant en oeuvre le procede 
de la presente invention, 

5 La Fig. 3 est un diagramme montrant les differentes etapes mises en oeuvre par un procede selon une premier 

mode de realisation de la presente invention, et 

La Fig. 4 est un diagramme montrant les differentes etapes mises en oeuvre par un procede selon une second 
mode de realisation de la presente invention. 

w [0031 J La presente invention est basee sur le fait que la matrice M (on rappelle que M = A+A dans le cas ZF ou M 
- A+A +• o 2 t dans le cas MMSE) representative du comportement d'un systeme de transmission CDMA est de type 
Toeptitz par blocs de dimensions [KxK]. 

[0032] En effet, si Ton regroupe les elements de la matrice Men blocs de dimensions KxK, les blocs qui appartiennent 
a une meme ligne paralleie a la diagonale allant du premier bloc ou dernier bloc de ia matrice M sont egaux entre eux. 
|| en est de meme des blocs qui appartiennent a ladite diagonale. 

[0033] On va considered en relation avec la Fig. 2 T une matrice U resultant de la decomposition de Cholesky d'une 
matrice M (U * C(M) « soit 0(A*A), soit C(A+A + <s% 

[0034] Comrne on peut le constater sur la Fig. 2, la matrice U est une matrice triangulaire superieure. Ainsi, chaque 
partie d'une ligne de blocs situee a gauche de la diagonale est constitute de blocs egaux a des matrices carres nulles 
20 de dimensions [K,K\. Quant a chaque partie d'une ligne de blocs d'indice i situee k la droite de la diagonale, diagonale 
incluse, elle est constitute d'au moins un bloc U jp de dimensions [K,K\ (p etant I'indice des colonnes de blocs) et au 
plus de P+ 1 blocs adjacents i de dimension [K,K] formant un ensemble de blocs, note par la suite E } . Les blocs qui 
sont situes au-dela du (P+lf dme , s'ils existent, sont des matrices nulles. On peut done constater a la Fig, 2 que : 
Ej - [U i0 U it1 ... U /(P / lorsque i<N-P+1 
25 E^IU] Q .!. U tN J\orsque /> N-P+1 et 

E N = U 

[0035] On rappelle que P etant la partie entiere de (1+(L-1)/Q). 

[0036] On peut noter egalement que I'element de la premiere ligne d'un ensemble de blocs E f peut egalement etre 
considere comme etant ie i,K dme fitment de la IK*™ 0 ligne de la matrice U. 
3o [0037] On notera qu'on utilise ici une notation ou I'indice de colonne de blocs p des blocs U lp est nul pour la diagonale 
de la matrice U. II s'ensuit que, dans la suite de la presente description, on dira qu'un element a n m (nft me ligne, rrf dme 
colonne) d'un bloc U jp d'un ensemble E f correspond a un element a n m (meme indices de ligne et de colonne) d'un 
bloc U Hp d'un ensemble E^. 

[0038] Tous les blocs sont differents les uns des autres. Cependant, J Rissanen a montre dans un article paru dans 
35 Mathematics of Computation, vol 27, n° 121 de janvier 1973, pages 147-157 qu'une matrice dont la dimension tend 
vers I'infini, qui est hermitienne, definie positive et de type Toeplitz en bloc et qui comprend un nombre fini de blocs 
non nul a pour matrice resultant de la decomposition de Cholesky une matrice U dont les ensembles E de blocs 
presentent des proprietes particulieres. 

[0039] En utilisant la notation de la Fig. 2, une de ces proprietes est que plus I'indice de lignes de blocs i est grand 
40 moins les blocs U ip (p compris entre 0 et P) sont differents des blocs U hi pi et ce, jusqu'fc tendre asymptotiquement 
vers I'egalite quand I'indice N tend vers I'infini. 

[0040] La presente invention est essentiellement basee sur cette propriete. La presente invention est egalement 
applicable du fait que la matrice a traiter dans le cas d'un reseau de telecommunication tel que celui qui est represents 
a la Fig. 1 est tres grande (par exemple dans le cas de I'UMTS-TDD ou A/vaut entre 61 et 976 et /Centre 1 et 16). 
45 [0041] La presente invention propose ainsi de calculer la matrice U de maniere iterative. A chaque iteration., un 
ensemble de blocs E, est calcule en appliquant ia methode de Cholesky et est ensuite ecrite dans la matrice U a sa 
place normale. 

[0042] Puis, sur la base de la valeur prise par un critere predetermine caracterisant revolution du processus decalcul 
de matrice U, il est decide soit du lancement de Titration sutvante i soit d'une duplication avec decalage a droite d'un 
50 bloc de la ligne d'indice / sur les autres lignes d'indice j avecj= M a N. Plus precisement, on considere une valeur 
prise par ce critere et, si elle est inferieure a une valeur seuil, la duplication avec decalage a droite d'un bloc est mise 
en oeuvre, sinon e'est iteration suivante qui est mise en oeuvre. 

[0043] Comme on le comprendra, lors de la duplication et du decalage a droite d'un bloc, il peut etre necessaire de 
tronquer I'ensemble a dupliquer lorsque le nombre de blocs a ecrire est inferieur au nombre de blocs de la ligne de 
55 blocs precedente. On notera que la valeur seuil pourrait varier dynamiquement en fonction des ressources disponibles 
au moment du traitement. 

[0044] Selon un premier mode de realisation, pour caracteriser revolution du processus de calcul de la matrice U, 
la valeur de ce critere depend des valeurs prises par au moins un element a nm de Pensemble E f de (a ligne de blocs 
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courante et le ou les Elements a nm correspondents de I'ensemble £ M de la ligne de blocs pr^cedente. 
[0045] Par exemple, ce critere a pour valeur la valeur absolue de ia difference des valeurs prises respectivement 
par, d'une part, ('element a n m de I'ensemble E } de la ligne de blocs courante choisi pour avoir la plus grande valeur 
absoiue parmi tous ses elements et pan d'autre part s Telement a n m correspondant de I'ensemble E H de la ligne de 
blocs precedents 

[0046] II peut egalement prendre pour valeur la valeur maximale des valeurs absolues des differences entre les 
valeurs prises respectivement par les Elements a n m de I'ensemble E/de la ligne courante et par les elements corres- 
pondants de I'ensemble EM de la ligne precedente. 

[0047] Dans ces deux cas, la valeur seuil correspond alors k un niveau d'erreurs acceptable. Le choix de sa valeur 
permet d'obtenir une detection conjointe qui est certes degradee (a un niveau correspondant k cette valeur) mais qui 
presente une complexite de calcul inferieure. 

[0048] On a repr&sente a la Fig. 3, un diagramme montrant le d^roulement du procedS de la presente invention 
selon ce premier mode de realisation. 

[0049] L'etape 1 est une etape ^initialisation pendant laqueile les differentes variables du processus sont determi- 
ners. II en estainsi notamment de la valeurs seuil S min et d'une variable /designant le numero d'ordre de iteration 
consideree qui est initialise ici a 1 . 

[0050] C'est a I'&ape 2 que sont calcuiees les valeurs de tous les elements de I'ensemble de blocs E v La ligne 
corrcspondante est memorisee dans la matrice U a I'endroit considere k une etape 3. 

[0051] A l'etape 4, sont sortis les elements de I'ensemble de blocs £ M qui vont servir a ia determination de la suite 
a donner au processus d'iteration. 

[0052] A l'etape 5, le numero cfordre de iteration i est increments, 

[0053] A l'etape 6, la valeur du critere S est calculee conformement k ce qui est decrrt ci-dessus. 
[0054] A l'etape 7, cette valeur S est comparee a la valeur seuil S min consideree d6terminee k l'etape 1 (ou even- 
tuellement recalculee depuis en fonction des conditions de ressources disponibles) et la decision quant a la suite du 
processus du calcul de la matrice Uest prise. Soit l'etape 2, est de nouveau mise en oeuvre pour le calcul de la ligne 
suivante. soit l'etape 8 de duplication et decalage a droite est mise en oeuvre. Pour ce faire, une boucle est prevue 
qui est interrompue lorsqu'& l'etape 9, il a 6te constate que la derniere ligne de ia matrice Ua 6te calcutee* 
[0055] Le resultat du precede de la ptesente invention est visible a la Fig. 2 ou ; a partir de la f dmo ligne, les lignes 
sont egales a cette ligne dupliquee puis deeaiee a chaque fois a droite d'un bloc, puis eventuellement tronquee lorsque 
la demi-ligne consideree comporte moins de P blocs. 

[0056] Selon un second mode de realisation, ce critere a pour valeur le nombre d'iterations realises. La valeur seuil 
correspond alors au nombre maximal d'iterations acceptables sans que soit d6pass£e la capacite de calcul maximale 
allouee au processus de detection conjointe. 

[0057] On a represents a la Fig. 4 : le diagramme correspondant a ce mode de realisation ou les memes etapes 
portent la meme reference et ou les etapes 4 et 6 ont disparu. L'etape T nouvelle effectue une comparaison entre le 
nombre i courant d'iterations et le nombre i max d'iterations acceptable tel que determine a l'etape 1 ou depuis selon 
les conditions de ressources disponibles. 

[0058] On notera que les premier el second modes de realisation ne s'excluent pas mutuellement ou pourraient etre 
mis en oeuvre conjointement. 

[0059] On notera encore que la description qui vient d'etre faite Pa 6te en considerant une decomposition de Cholesky 
d'ou resultait une matrice triangulaire superieure Uei que ia prSsente invention s'applique egalement a une decom- 
position de Cholesky conduisant a une matrice triangulaire inferieure L Le procSde comporterait alors les memes 
etapes que celui qui vient d'etre decrit. 



Revendications 

1 . Procede de detection conjointe destine a etre mis en oeuvre dans un rScepteur d J un systeme de transmission, te! 
qu'un systeme de transmission k acces multiple a repartition de codes, ladite detection consistant en une etape 
de calcul des signaux qui ont ete emis par lesdits emetteurs dudit systeme mis sous ia forme d'un vecteur xet ce, 
a partir de I'ensemble des signaux re<?us des Smetteurs dudit systeme mis sous la forme d'un vecteur y, de I'en- 
semble des caracteristiques dudit systeme mis sous forme d'une matrice A et Sventuellement de I'ensemble des 
caracteristiques de bruit mis sous forme d'une matrice n, ladite etape de calcul consistant en une etape ^inversion 
d'une matrice M fonction de la matrice systeme A, une etape d'application de I'ensemble des signaux re$us y sur 
une matrice resultant du produit de la matrice inverse Itf 1 par le transpose conjuguS de la matrice A, ladite etape 
d'inversion de la matrice M consistant en sa decomposition de Cholesky en une matrice triangulaire Uou L, ca- 
racterise en ce qu'il consiste a determiner ladite matrice triangulaire de maniere iterative, chaque iteration consis- 
tant d'abord au calcul d'un ensemble de blocs appartenant a une ligne de blocs de ladite matrice puis, en fonction 
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de la valeur prise par un critere predetermine caracterisant revolution du processus de calcul de ladite matrice 
decomposee comparee a une valeur seuil, soit au passage a Iteration suivante soil a la duplication de la ligne 
courante sur toutes les autres lignes restantes avec decalage a droite d'un bloc et, eventuellement, tronquage des 
blocs en excedant. 

Procede selon la revendication 1 . caracterise en ce que ladite valeur seuil est determinee dynamiquement en 
fonction des ressources disponibies au moment du calcul. 

Procede selon la revendication 1 ou 2 : caracterise en ce que ledit critere est tel que sa valeur depend des valeurs 
respectivement prises par au moins un element de I'ensemble de blocs de la lighe courante et le ou les elements 
correspondants de i'ensemble de blocs de la ligne calculee k Iteration precedente. 

Procede selon la revendication 3, caracterise en ce que ledit critere est tel que sa valeur est egale a la valeur 
absolue de la difference des valeurs prises respectivement par Telement de I'ensemble de blocs de la ligne courante 
dont la valeur absoiue est la plus grande et par Telement correspondant de I'ensemble de blocs de la ligne pre- 
cedente. 

Procede selon la revendication 3, caracterise en ce que ledit critere est tel que sa valeur est egale 'd la valeur 
maximal© des valeurs absolues des differences entre les valeurs prises respectivement par tous les elements de 
I'ensemble de blocs de la ligne courante et par les elements correspondants de I'ensemble de blocs de la ligne 
precedente. 

Procede selon une revendlcations precedentes, caracterise en ce que ledit critere a pour valeur le nombre d'ite- 
rations, la valeur seuil correspondant aiors a un nombre maximal d'iterations. 

Recepteur d'un systeme de transmission, tel qu'un systeme de transmission a acces multiple a repartition de 
codes, prevu pour recevoir et detecter les signaux qui ont ete emis par des emetteurs dudit systeme de transmission 
afin de les recouvrir, caracterise en ce que, pour detecter iesdits signaux, il est prevu pour mettre en oeuvre un 
procede de detection selon une des revendications 1 a 6. 
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The present invention relates to a joint detection method 
intended to be applied in a receiver of a transmission system, 
such as a transmission system of the code division type CDMA 
5 (Code Division Multiple Access) . 

A large number of telecommunications systems being 
developed apply CDMA techniques . These techniques use code 
sharing of frequency resources, i.e., several users use the 
same frequency band at the same time, and are differentiated 

10 by different spreading codes which they use. A user generally 
distinguishes the signal which is addressed to him, in the 
total received signal, by comparing the latter with his/her 
spreading code. 

However, the quality of the received signal is reduced by 

15 residual interference from other users called MAI (Multiple 
Access Interference) . Various techniques including joint 
detection, more particularly used in the case of synchronous 
spreading codes, are directed to using the knowledge of 
spreading codes of the other users in order to reduce this 

20 residual interference from the other users. 

A model of the upstream link of a CDMA type transmission 
system which uses a joint detection method, is illustrated in 
Fig. l.This system is provided for receiving K users which 
provided with a transmitter 10 k (k = 1 to K) . Each transmitter 

25 10 k is provided for transmitting data d k of a user k by bursts 
of N data symbols, these data being then noted as d k = fx 1 *, 
wherein x"* then represents the nth symbol transmitted by 
the user k. In each transmitter 10 k , the data d k of a user k 
are submitted to a spreading code c k with a length of Q bits 

30 allocated to the relevant user k before being transmitted over 
a channel 20 k . Each channel 20 k is characterized by an impulse 
response varying over time noted as h k and a frequency band 
width noted as B c = 1/T C . The data are then transmitted over a 
channel 20 k . For estimating the impulse response h k of each 



2 



10 



15 



20 



channel 20 k , L samples spaced out from each other by a time 
interval equal to T c are taken into account. 

Klein et al . , have been able to show that the signal y 
received by the receiver may be represented by a vector y of 
dimension (N X Q + L - 1) , the components of which are 
representative of the transmitted signals by each of the K 
users, possibly added to Gaussian noise represented by a 
vector n also of dimension (N x Q + L-l) . So we may write: 

y = Ax + n 

wherein 

X I X . X m • • * « X » * " * ■ X , X m m m • X I 

is the signal transmitted before spreading, 
is the signal received by the receiver, 



is the additive Gaussian white noise, and 

A is a matrix of dimensions (N. Q + L-l) X (K.N) 
Each element of matrix A, noted as A±j, is determined 
25 according to the estimation of the impulse response of each of 

the channels of order k: 



30 



and according to each of the spreading codes of order k: 



3 



c'=[c 1 *,c*,..,c*]avect = l,..., J fi: 



If one sets 



i=Q.(n-1)+ w 



10 



15 



and 



j = k+K(n-1) 

A CDMA type transmission system may be modeled in such a 
way that each element Aij may be written as: 

A n = b k for *=1..Kin=1...Wand iv=1..Q*H 

Ajj=Q else 

with 

b k =c k *h k ..,b* +L _ t ] for k = \..X 



In order to estimate the data transmitted by each user k, 
one may proceed in receiver 30, with a joint detection such as 
the one described by A. Klein, G.K. Kaleh and P.W. Baier in an 
article entitled "Equalizers for multi-user detection in code 
20 division multiple access mobile radio systems" published in 
IEEE of 1994 . 

These known joint detection methods essentially consist 
of solving the matrix equation given above by determining the 
vector of symbols x by knowing the vector of the detected 
25 signal y, the system matrix A and the noise matrix n. 

It will be noted that matrices y, A and n, are moreover 
determined by estimation methods which are known to one 



4 



skilled in the art and they are therefore not described in 
more detail in the present specification. 

It will still be noted that, notably from the fact that 
the impulse response of each of the channels 20 k varies over 
5 time, it is at each instant that the elements of matrix A 
change values . 

Determination of the vector of symbols x± by the method 
expounded above generally uses one of the two following 
criteria which are zero forcing or ZF on the one hand, and on 
10 the other hand the minimum mean square error (MMSE) . Both of 
these criteria respectively lead to the following solutions: 

x = (A+AY 1 A+ v 

{ ' J for the zero forcing criterion 

and 

x = (A+A+fc^iy^A+v 

K K ' 1 for the minimum mean square error MMSE 

15 criterion (J is the identity matrix, d 2 being the variance of 
the elements of the noise matrix n) . 

In the following description, the exponent + applied to a 
matrix means that it is its conjugate transpose. 

These equations may also be written as: 

20 

x=M'Vy 



As it may be seen, solving this equation therefore 
requires inverting a matrix noted as M, which is either equal 

25 to the product A*. A of the conjugate transpose system matrix A + 
by the system matrix A, or equal to the product A* .A+o 2 1 of the 
conjugate transpose system matrix A + by the system matrix A to 
which is added the identity matrix I multiplied by the 
variance a 2 of the noise. 

30 When it is hermitian, defined and positive, (as this is 

always the case with the modeling of a CDMA system given 
above) , a matrix M may be decomposed, by means of an 



5 



algorithm, the so-called Cholesky algorithm, into an upper 
triangular matrix U (or lower triangular matrix L according to 
the chosen notation) such that matrix M is equal to the 
product of the conjugate transpose upper triangular matrix if 
5 by matrix U: 

M=lf. U 

or depending on the case, to the product of the lower 
10 triangular matrix L by the conjugate transpose matrix L + ': 

M=LZ_ + 

now, inversion of matrix M is largely simplified as the 
15 inverse matrix XT 1 is either equal to the product of the 
inverse matrix IT 1 by the conjugate transpose of this matrix 

rrrV: 



M' 1 = U 1 (If 1 )* 

20 

either to the product of the conjugate inverse matrix (L~ 
1 ) + by matrix IT 1 : 



M' 1 = (L' 1 ) + L' 1 

25 

Because of the few operations required for inverting U, 
inversion of M by passing through the Cholesky decomposition 
requires much less operations than inversion of M by direct 
calculation . 

30 This method with Cholesky decomposition allows the 

complexity of this calculation to be reduced by a factor of 
about 2 (complexity of T 3 /6, T being the dimension of the 
matrix) as compared with conventional direct inversion 
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techniques. However, the amount of calculations still remains 
very significant. 

The object of the present invention is to provide a link 
detection method which uses the method with the Cholesky 
5 decomposition but for which the amount of calculations is 
reduced as compared with methods from the prior art. 

For this purpose, a link detection method according to 
the present invention is of the type which consists of a step 
for calculating the signals which have been transmitted by 

10 said transmitters of said system formed as a vector x and this 
from the set of signals received from the transmitters of said 
system, formed as a vector y, from the set of characteristics 
of said system, formed as a matrix A and possibly from the set 
of noise characteristics formed as a matrix n, said 

15 calculation step consisting of a step for inverting a matrix M 
depending on the system matrix A, a step for applying the set 
of received signals y on a matrix resulting from the product 
of the inverse matrix AT 1 by the conjugate transpose of matrix 
A, said step for inverting matrix M consisting in its Cholesky 

20 decomposition into a triangular matrix U or L. 

This method is characterized in that it consists of 
determining said triangular matrix iteratively, each iteration 
first consisting of calculating a set of blocks belonging to a 
line of blocks of said matrix and then, depending on the value 

25 assumed by a predetermined criterion characterizing the 
development of the process for calculating said decomposed 
matrix as compared with a threshold value, either by passing 
to the next iteration or by duplicating the current line on 
all the other remaining lines with right-shifting of a block 

30 and possibly truncation of excess blocks. 

According to another feature of the invention, said 
threshold value is dynamically determined depending on the 
available resources at the moment of the calculation. 

According to another feature of the invention, said 

35 criterion is such that its value depends on values 
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respectively assumed by at least one element of the set of 
blocks of the current line, and the corresponding element (s) 
from the set of blocks of the line calculated in the previous 
iteration. 

5 According to another feature of the invention, said 

criterion is such that its value is equal to the absolute 
value of the difference of values assumed by the element of 
the set of blocks of the current line, for which the absolute 
value is the largest, and by the corresponding element of the 
10 set of blocks from the previous line, respectively. 

According to another feature of the invention, said 
criterion is such that its value is equal to the maximum value 
of the absolute values of the differences between the values 
assumed by all the elements of the set of blocks of the 
15 current line and by the corresponding elements of the set of 
blocks from the previous line, respectively. 

According to another feature of the invention, said 
criterion has the value of the number of iterations, the 
threshold value then corresponding to a maximum number of 
20 iterations. 

The present invention also relates to a receiver which is 
provided for applying a detection method such as the one which 
has just been described. 

The aforementioned features of the invention, as well as 
25 other features, will become more apparent upon reading the 
following description of exemplary embodiments, said 
description being made in connection with the appended 
drawings, wherein: 

Fig. 1 is a block diagram of a CDMA system which is 
30 provided for applying a joint detection method according to 
the present invention, 

Fig. 2 is a diagram of a matrix resulting from a Cholesky 
decomposition applying the method of the present invention. 
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Fig. 3 is a flow chart showing the different steps 
applied by a method according to a first embodiment of the 
present invention, and 

Fig. 4 is a flow chart showing the different steps 
5 applied by a method according to a second embodiment of the 
present invention . 

The present invention is based on the fact that the 
matrix M (let us recall that M = A + A in the ZF case, or M = A + A 

10 + d*I in the MMSE case) representative of the behavior of a 
CDMA transmission system is of the Toeplitz type per blocks of 
dimensions [K X K] . 

Indeed, if the elements of matrix M are grouped into 
blocks of dimensions KX K, the blocks which belong to a same 

15 line parallel to the diagonal ranging from the first block to 
the last block of the matrix are equal to each other. The same 
applies to blocks which belong to said diagonal. 

In conjunction with Fig. 2, a matrix U resulting from the 
Cholesky decomposition of a matrix M ( U = C(M) = either 

20 C(A+A), or C (A + A) + o 2 1) will be considered. 

As it may be seen in Fig. 2, matrix U is an upper 
triangular matrix. Thus, each portion of a line of blocks 
located on the left of the diagonal consists of blocks equal 
to zero square matrices with dimensions [K,K] . As for each 

25 portion of a line of blocks with index i located on the right 
of the diagonal, including the diagonal, it consists of at 
least one block U± /P of dimensions [K,K] (p being the index of 
the block columns) and at the most of P+l adjacent blocks i of 
dimensions [K,K] forming a set of blocks, noted as Ei 

30 subsequently. The blocks which are located beyond the (P+l) th 
one, if they exist, are zero matrices. So, it may be seen in 
Fig. 2 that: 
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E i = f u i,o u i, 1 - u i tP J when / <W-P-f / 
E, = [U, 0 ... U iN .ij when i> N-P+1 and 
E N =U 

Let us recall that P is the integer portion of (1+(L- 
D/Q) . 

It may also be noted that the element of the first line 
5 of a set of blocks E± may also be considered as the i.K th 
element of the i. inline of matrix [7. 

It will be noted that a notation is used here where the 
column index of the p blocks from the U± fP blocks is zero for 
the diagonal of matrix U. It follows that, in the following of 
10 the present description, it will be stated that an element a n , m 
(n th line, m th column) of a block Ui fP of a set E± corresponds 
to an element a n , m (same line and column indexes) of a block U±- 
ltP of a set Bi-i. 

All the blocks are different from each other. However, J. 
15 Rissanen shows in an article published in Mathematics of 
Computation, Vol. 27, No. 121, January 1973, p 147-157, that a 
matrix with a dimension tending to infinity, which is 
hermitian, defined positive and of the Toeplitz block type and 
which comprises a finite number of non-zero blocks has as 
20 resulting matrix from the Cholesky decomposition, a matrix U 
for which the sets E of blocks have particular properties . 

By using the notation of Fig. 2, one of these properties 
is that the larger the index of the lines of blocks i, the 
lesser the blocks U± fP (p between 0 and P) are different from 
25 blocks Ui- lfP > and this until they asymptotically tend to 
equality when index N tends to infinity. 

The present invention is essentially based on this 
property. The present invention is also applicable because the 
matrix to be processed in the case of a telecommunications 
30 network such as the one illustrated in Fig. 1 is very large 
(for example, in the case of UMTS -TDD where N is between 61 et 
96 and K between 1 and 16) . 
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The present invention thus suggests that the matrix U be 
calculated iteratively. At each iteration, a set of blocks Ei 
is calculated by applying the Cholesky method and it is then 
written into the matrix U at its normal place. 
5 Next, on the basis of the value assumed by a 

predetermined criterion characterizing the development of the 
process for calculating matrix U, a decision is made whether 
to initiate the following iteration or to duplicate a block 
with right-shifting from the line of index i onto the other 
10 lines of index j with j = i+1 to N. More specifically, a value 
assumed by this criterion is considered and if it is less than 
a threshold value, duplication of a block with right-shifting 
is applied otherwise it is the following iteration which is 
applied . 

15 As it will be understood, upon duplicating and right- 

shifting a block, it may be necessary to truncate the set to 
be duplicated as the number of blocks to be written is less 
than the number of blocks of the previous line of blocs. It 
will be noted that the threshold value may dynamically vary 

20 according to the available resources at the moment of the 
processing . 

According to a first embodiment, to characterize the 
development of the process for calculating matrix U, the value 
of this criterion depends on values assumed by at least one 

25 element a n/m of the set Ei of the current line of blocks and 
the corresponding element (s) a n , m of the set E±-i of the 
previous line of blocks. 

For example, this criterion has the value of the absolute 
value of the difference of values assumed by element a n , m of 

30 the set Ei of the current line of blocks selected in order to 
have the largest absolute value among all its elements on the 
one hand, and on the other hand by the corresponding element 
an,m of the set Ei-i of the previous line of blocks, 
respectively. 
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It may also have the value of the maximum value of the 
absolute values of the differences between the values assumed 
by the element a n , m of the set E± of the current line and by 
the corresponding elements of the set Ei- ± of the previous 
5 line, respectively. 

In both of these cases, the threshold value then 
corresponds to an acceptable error level. By selecting its 
value, it is possible to obtain a joint detection which is 
certainly degraded (at a level corresponding to this value) 
10 but which has a lesser complexity of calculation. 

A flow chart showing the development of the method of the 
present invention according to this first embodiment is 
illustrated in Fig. 3. 

Step 1 is an initialization step during which the 
15 different variables of the process are determined. This is 
notably the case of the threshold values S m in and of the values 
of a variable i designating the order number of the relevant 
iteration which is initialized to 1 here. 

It is in step 2 that the values of all the elements of 
20 the set of blocks Ei are calculated. The corresponding line is 
stored in matrix U at the location considered in a step 3 . 

In step 4, the elements of the set of blocks Ei-± are 
outputted, which will be used for determining the follow-up as 
to the iteration process. 
25 In step 5, the order number of iteration i is 

incremented. 

In step 6, the value of criterion S is calculated 
according to what is described above. 

In step 7, this value S is compared with the relevant 
30 threshold value S m ± n as determined in step 1 (or possibly as 
re-calculated since then, depending on the conditions of 
available resources) and the decision is made as to the 
follow-up of the process for calculating the U matrix. Either 
step 2 is again applied for calculating the next line, or step 
35 8 for duplication and right-shifting is applied. To do this, a 
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loop is provided which is interrupted when in step 9, it has 

been seen that the last line of matrix U has been calculated. 

The result of the method of the present invention is 

visible in Fig. 2 where, from the j th line, the lines are equal 
5 to this line which is duplicated then right-shifted every time 

by one block, and then possibly truncated when the relevant 

half -line includes less than P blocks. 

According to a second embodiment, this criterion has the 

value of the number of achieved iterations. The threshold 
10 value 'then corresponds to the maximum number of acceptable 

iterations without exceeding the maximum calculation capacity 

allocated to the joint detection process. 

The flow chart corresponding to this embodiment where the 

same steps bear the same reference and where steps 4 and 6 
15 have disappeared, is illustrated in Fig. 4. The new step 7' 

carries out a comparison between the current number i of 

iterations and the acceptable number of iterations, as 

determined in step 1 or since then, according to the 

conditions of available resources. 
20 It will be noted that the first and second embodiments 

are not mutually exclusive or may be applied together. 

It will further be noted that the description which has 

been made, was made by considering a Cholesky decomposition 

from which an upper triangular matrix U resulted and that the 
25 present invention also applies to a Cholesky decomposition 

leading to a lower triangular matrix L. The method would then 

include the same steps as the one which has been described. 
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